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In this paper, transformation optics is presented together with a generalization of invisibility cloaking: instead of an empty 
region of space, an inhomogeneous structure is transformed via Pendry's map in order to give, to any object hidden in the central 
hole of the cloak, a completely arbitrary appearance. Other illusion devices based on superlenses considered from the point of view 
of transformation optics are also discussed. 
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I. Introduction 

In 2006, it was suggested by Pendry et al. [1] that an object 
surrounded by a coating consisting of an exotic material could 
become invisible to electromagnetic waves. This device was 
named "invisibility cloak" in reference to Harry Potter, the 
popular character of J.K. Rowling. Beside his famous cloak, 
the little wizard has other spells to go unnoticed. Among 
the most spectacular is the "polyjuice potion" that is able 
to turn somebody into anybody else's appearance [2]. In this 
paper, we do not present a potion but rather an optical device 
able to accomplish the same task, i.e. to give an arbitrary 
optical response chosen in advance to any other object placed 
inside the device. In fact, the principle is here very similar 
to the design of Pendry's invisibility cloak but, instead of 
geometrically transforming an empty domain, we transform 
a region containing the object to be imitated, thus leading to 
a generalization of cloaking. 

II. Transformation Optics 

In recent years, transformation optics has become a very 
active new field. It has been popularized through the idea of 
J.B. Pendry that an invisibility cloak can be designed by trans- 
forming space and considering the corresponding equivalent 
material properties [1], [3]. Indeed, it is a deep property of 
Maxwell's equations that they are purely topological (when 
written in the proper formalism [4]) and that all the metric 
aspects can be encapsulated in the electromagnetic mate- 
rial properties. A direct consequence is that any continuous 
transformation of space can be encoded in an equivalent 
permittivitty and permeability. Extending this principle beyond 
continuous transformations allows to design exotic optical 
devices such as the invisibility cloak. 

Exterior calculus is the most natural formalism to write 
Maxwell's equation [5], [6] so that they have the following 
form (in the harmonic case with a pulsation to and complex 
valued fields): 

dH = J — iujT) 
dE = +iu>B 
dT> = p 
dB = 



(1) 



where d is the exterior derivative (d plays the role of curl in 
the two first equations and of div in the last two ones, see 
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Appendix), the 1 -forms E, H are the electric and magnetic 
fields respectively, the 2-forms D, B, and J are the electric 
flux density or displacement, the magnetic flux density or 
induction, and the electrical current density respectively, and 
the 3-form p is the electric charge density. The only operator 
involved is the exterior derivative that is completely indepen- 
dent from the metric. 

The metric is involved in the Hodge star operator * (see 
Appendix) that is necessary to introduce the constitutive laws 
of materials (including a void, where, D = £o * E and B = 
/io * H). It can also be argued that it is in fact these very 
electromagnetic properties of space that determine the metric 
[7]. This formalism has proved to be very useful in the context 
of the numerical solution of Maxwell's equations [4], [8]. In 
this case, in has been shown that the topological structure of 
the equation can be preserved at the discrete level (for instance 
in Yee's FDTD algorithm or using Whitney discrete forms as 
finite elements) while the whole process of approximation is 
concentrated in the design of the discrete Hodge operator [9]. 

A. Change of Coordinates in Maxwell's Equations 

In the exterior calculus formalism, the only task associated 
to changing a coordinate system is to determine an explicit 
expression for the Hodge star operator [4], [10]. A very useful 
point of view is to consider weak formulations where integrals 
of volume forms (3-forms) are built with scalar products of 
forms, i.e., exterior products together with the Hodge operator 
acting on one of the factors. 

For instance, the wave equation for the electric field (in 
homogeneous media): 



d(n 1 * dE) - oj 2 e * E 



0. 



(2) 



has the following weak formulation: find E e i7(curl, O) 
such that 



J n pT x * dE A dE dx - cv 2 J Q e * E A E dx = , 

VE' e ff (curl,n) 



(3) 



where A is the exterior product (see Appendix). 

We can use the fact that we know how to write this ex- 
pression components by components in a Cartesian coordinate 
system and that we also know how to transform the derivative 
and the multiple integrals to determine the action of the Hodge 
operator in other coordinate systems. 
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Considering a map from the coordinate system {u, v, w} 
to the coordinate system {x, y, z} given by the functions 
x(u, v,w), y(u,v,w), and z(u, v, w), the transformation of 
the differentials is given by: 



with 



\dy=^du+^dv + ^d W 
{dz=^du+^dv + ^dw. 



(4) 



Given a p-form expressed in the {x, y, z} coordinate system, 
it suffices to replace the dx,dy, dz by the corresponding 1- 
forms involving du, dv, dw in the basis exterior monomials to 
obtain the expression of the form in the new coordinate system. 
Note that the form travels naturally counter to the current with 
respect to the map and this is why this transportation of the 
forms from x, y, z to u, v, w is called a pull-back. 

This operation can be defined not only between two coor- 
dinate systems but also between two different manifolds even 
if they do not have the same dimensions. 

Consider two manifolds (or more simply, two open domains 
of R m and E" respectively) N and M and a (regular) map 
(p from N to M such that (for simplicity) ip(N) = M. 
The example above shows that it is very easy to express 
the differentials of the coordinates on M in terms of the 
differentials of the coordinates on N and therefore to find the 
image on N of a 1-form on M given by the dual map tp*, from 
M to N, also called, as indicated above, the pull-back. In fact, 
any covariant object such as a p-form or a metric can be pulled 
back by translating the differentials on M into differentials 
on N. Defined in this way, the operation commutes of course 
with the exterior and tensor products but also with the exterior 
derivative and the Hodge star (defined with the pulled-back 
metric) [4]. 

As for contravariant objects such as vector fields, they travel 
forward just like the geometrical domains. Given a vector v at 
a point p on N, it suffices to choose a curve 7 going through 
the point and such that the vector is the tangent vector to the 
curve at this point, to take the image of the curve ^(7) on M 
and the vector tangent to this curve at the point <p(p) as the 
image of v. Defined in this way, the map for vectors from N to 
M, denoted by </9*(v) or dip(v), is called the differential of <p 
or the push-forward and it can be extended to any contravariant 
object. 

Another fundamental property of the pull-back is its com- 
mutativity with integration in the sense that, for any form a 
that is integrable on a subset <p(d) of M, which is the image 
of a subset f2 of N, one has: 



/ a = / ip*(a). 



(5) 



All the information for the pull-back is therefore contained 
in the Jacobian matrix J (or maybe we should say matrix field 
since it depends on the point in space considered) in terms of 
which Eq. (4) can be written: 



(6) 



dx \ 




dy 


I = J ( dv 


dz ) 


\ dw 



3(u, v, w) 



d(x,y,z) 
d(u, v, w) 



dx 


dx 


dx 


du 
dy 




dw 
dy 






dw 


<& 


I 


dz 


du 


dv 


dw 



Using matrix notation, the detailed computation of the 
relation between the coefficients of a 1-form a in {x,y,z} 
and {u, v, w} coordinates is performed as follows: 

(dx 
dy 
dz 

du 

= (a x a y a z ) J ( dv 
dw 

du 

a u du + a v dv + a w dw = (a u a v a w ) | dv 

dw 



and the following relation is obtained: 

(a x a y a z ) J = (a u a v a w ). 



(7) 



Now the contributions to weak form integrals like Eq. (3) may 
have the following form: 



/ 

Jn 



a A *a\ 



where a and a' are 1 -forms (that can be obtained as gradients 
of a scalar field, although it really does not matter here). 
The question is: how to deal with the Hodge operator? A 
direct attack would be to pull back the metric and use the 
explicit expression of the operator but it is faster here to take 
advantage of the simple form of the scalar product in Cartesian 
coordinates that reduces to the dot product. Again using matrix 
notation (where 3 T is the inverse of J T ): 

a A *a' = (a x a y a z )(a' x a' y a' z ) T dx A dy A dz 

= (a u a v a M ,)J _1 {(a' u a' v a' w ) J" 1 ] 7 ^ A dy A dz 

= (a u a v a w )J~ 1 J~ T (a' u a' v a' w ) T det(3)du A dv A dw. 

(8) 

The first line is the definition of the scalar product of 1-forms 
equated to the scalar product in Cartesian coordinates. 

The fact that the transformation of 3-forms dx Ady Adz = 
det(J)c?u A dv A dw only involves the Jacobian, i.e. the 
determinant of the Jacobian matrix, has been used here. Hence, 
the only difference from the case of Cartesian coordinates is 
that one of the (column) vectors has to be multiplied on the 
left by a symmetric matrix T _1 before performing the dot 
product, where T is given by: 

J T J 

~ dct(J) ' 



(9) 



It is now interesting to look at how a particular 2-form basis 
monomial transforms, for instance, 

,dx , dx , dx , n r dy , dy , dy , . 
dxAdy = \—du+—dv+—dw]A\—du+—dv+—dw] 
du dv dw du dv dw 

/ dx dy dx dy \ ^ ^ / dx dy dx dy\ ^ ^ 
\ du dv dv du J \ dv dw dw dv J 
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dw A du. 



(10) 



dx dy dx dy 
dw du du dw , 

The cofactors of J are now involved in the transformation. 
(These are the elements of J~ T dct(J).) 
Given a 2-form: 

f3 = fi x dy A dz + (3 y dz A dx + f3 z dx A dy 
= j3 u dv Adw + (3 v dw Adu + (3 w du A dv 

the following relation is obtained: 

J- T det(J) = (/3 u /3 v p w ) (11) 

and, considering the scalar product j3 A *f3' of two such 2- 
forms, it is straightforward to show that the matrix involved 
in the transformation of this scalar product is here T (instead 
of its inverse in the case of the scalar product of two 1-forms): 

P A */3' — [fi x /3 y p y (i' z ) T dx A dy A dz 

= ifiu PvMag^ [(/? « & ^ W 1 T det W du A (iw A dw 



dot(J) V^m 



3' ^ 

v l w 



) du A dv A dw. 



(12) 

Everything can now be summarised in the following recipe 
that takes into account implicitly the Hodge star: consider a 
3-form 7 to be integrated on a domain O in order to get J n 7 
to contribute to a weak form, then: 

• If the integrand involves only scalars (0-forms or 3-forms 
and it does not matter if the 3-forms are expressed as 
the divergence of a vector field) or if it is the exterior 
product of a 1-form and a 2-form (and it does not matter 
if they are respectively a gradient and a curl) looking 
superficially like a scalar product of vectors, only det(J) 
has to be introduced as a factor. 

• If the integrand is the scalar product of two 1-forms (and 
it does not matter if one or both 1-forms are expressed 
as the gradient of a scalar field), multiply on the left one 
of the column vectors of coefficients by the matrix T _1 . 

• If the integrand is the scalar product of two 2-forms (and 
it does not matter if one or both 2-forms are expressed 
as the curl of a vector field), multiply on the left one of 
the column vectors of coefficients by the matrix T. 

The expression obtained for if*(j) depending on variables 
u, v and w (coordinates x, y and z have been replaced by the 
functions x(u, v,w), y(u,v,w) and z(u, v, w) respectively) is 
integrated on ft to get the desired contribution to the volume 
integral of the weak formulation. 

It can also be interesting to consider a compound transfor- 
mation, i.e. the transformation of a transformation. Consider 
three systems of coordinates u i7 Xi and Xi (possibly on 
different manifolds) and the maps ip Xu '■ Ui -> Xi given by 
functions Xi(uj) and (p xX \ Xi — > Xi given by functions 
Xi(Xj). The composition map (p xX o ip Xu — fxu '■ — > is 
given by the functions: Xi(Xj(uk)). If 3 xX and J Xu are the 
Jacobian matrices of the maps <p xX and ip Xu respectively, the 
Jacobian matrix J xu of the composition map ip xu is simply 
the product of the Jacobian matrices: 

^xu = xX 3 Xu . 

This rule naturally applies for an arbitrary number of maps. 

It is also worth noting that the matrix J T J is nothing but 
the metric tensor whose coefficients are expressed in the local 
coordinates. 



B. The Geometric Transformation - Equivalent Material Prin- 
ciple 

A very interesting interpretation of the preceding formulae 
is that the matrix T and its inverse can be viewed as tensorial 
characteristics of equivalent materials. 

By inspection of Eq. (3), it appears that /i -1 is present as 
a factor in the term involving the exterior derivatives, i.e. a 
scalar product of two 2-forms and that the T factor can be 
introduced by multiplying /i by T _1 (and therefore turning 
it in a tensor quantity). It appears also that e is present as a 
factor in the term involving directly the electric field, i.e. a 
scalar product of two 1-forms and that the T^ 1 factor can be 
introduced by multiplying e by T _1 (and therefore turning it 
also in a tensor quantity). 

Therefore, the only thing to do in the transformed co- 
ordinates to compute the integrals of the weak form is to 
replace the materials (often homogeneous and isotropic) by 
equivalent ones that are inhomogeneous (their characteristics 
are no longer piecewise constant but merely depend on u, v, w 
coordinates) and anisotropic ones (tensorial nature) whose 
properties are given by 



eT 



and n' = nT^ 1 . 



(13) 



We note that there is no change in the impedance of 
the media since the permittivity and permeability undergo 
the same transformation. As for the vector analysis operator 
and product, everything works as if we were in Cartesian 
coordinates. 

In electromagnetism, changing a material can thus be 
viewed as changing metric properties and conversely a change 
of coordinates can be taken into account by introducing a 
fictitious equivalent material. For a general transformation, 
the equivalent material is inhomogeneous and anisotropic. It 
may be interesting in some cases to introduce non-orthogonal 
coordinate systems to facilitate the solution of particular 
problems, e.g. helicoidal geometries [11], [12]. 

It is straightforward to generalize the present rules to 
initially anisotropic material properties. It has also to be noted 
that they need not to be initially homogeneous. Therefore, 
the basic principle of transformation optics can be stated in a 
very general setting: For all our practical purposes, M and 
will be here the whole or parts of W 3 . Given a map 
(f from a space to a space M determining a geometric 
transformation (i.e. given a Cartesian coordinate system x on 
M and an arbitrary coordinate system x' on N, ip : N — >• M 
is described by x(x'), i.e. x given as function of x'), when 
one has an electromagnetic system described by the tensor 
fields e(x) for the dielectric permittivity and /x(x) for the 
magnetic permeability in the space M, if one replaces the 
initial material properties by equivalent material properties 
given by the following rule [10], [13], [11], [12]: 



gV) = J- 1 (x')£(x(x'))J" T (x')dct(J(x')), 
/x'(x') = J- 1 (x')M(x(x'))J" T (x')dct(J(x')), 



(14) 



one gets an equivalent problem on N. Here, an equivalent 
problem means that the solution of the new problem on N, 
i.e. electromagnetic quantities described as differential forms, 
are the pulled back of the solution [4] of the original problem 
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on M and that the same Maxwell's equations (i.e. as if we 
were in Cartesian coordinates or, more accurately, having the 
same form as (1) written with the exterior derivative) are still 
satisfied. 

The transformation rules between M and TV, for compo- 
nents of the fields that are 1-forms such as E and H, is 
given by Eq. (7) and, for the components of the fields that 
are 2-forms such as B, D, and J, is given by Eq. (11). It 
appears clearly that the invariant quantities, according to Eq. 
(5), are the global quantities built as integral of p-forms on 
p dimensional geometrical objects: the line integrals of the 
electric field and of the magnetic field along a curve, the fluxes 
across a surface of the electric displacement, the magnetic flux 
density, the current density, the Poynting vector... 

III. Generalized cloaking 

In this section, we present a generalization of cloaking able 
to arbitrarily transform the electromagnetic appearance of an 
object. The basic principle is to obtain the constitutive relations 
of the cloak by application of a space transformation to a non- 
empty region. Invisibility can be considered as a particular 
case that corresponds to choosing the empty space as the object 
to be faked. 

In the case of the cylindrical Pendry's map [1], [3], [14], 
described by the transformation of the 2D cross section, the 
plane R 2 minus a disk Di of radius i?i is mapped on the whole 
plane R 2 in such a way that a disk D 2 of radius R 2 > Ri, 
concentric with D\, is the image of the annulus D 2 \D\ by a 
radial transformation (see Fig. 1). In cylindrical coordinates, 
this transformation is given by: 

jr = (r' - R 1 )R 2 /(R 2 - Ri) for R x < r' < R 2 , 

\0 = 0>,z = z'. U; 

As for the outside of the disk D 2 , the map between the two 
copies of R 2 \D 2 is the identity map. 

The material properties given by rule (14) corresponding to 
this transformation provide an ideal invisibility cloak: outside 
D 2 , everything behaves as if we were in free space, including 
the propagation of electromagnetic waves across the cloak, and 
is completely independent of the content of D\. 

Now, rule (14) may be applied to D 2 containing objects with 
arbitrary electromagnetic properties so that a region cloaked 
by this device is still completely hidden but has the appearance 
of the objects originally in D 2 . We may call this optical effect 
masking [15] or "polyjuice" effect. 

IV. Numerical modeling 

Figs. 4 and 5 show the effect of masking on a scattering 
structure. On Fig. 4, a cylindrical TM wave emitted by 
a circular cylindrical antenna is scattered by a conducting 
triangular cylinder (the longest side of the cross section is 
1.62A and e r = 1 + 40i). The field map represents the 
longitudinal electric field E z (x, y) and the outer boundary of 
the cloak is shown to ease the comparison with the masked 
case. On Fig. 5, the same cylindrical TM wave is scattered by a 
masked triangular cylinder (but the scattering object inside the 
cloak may be arbitrary as far as it is small enough to fit inside 




Initial (Cartesian coordinates) New 



Fig. 1. Pendry's map of an annulus to a disk used to determine the material 
properties of an invisibility cloak via the equivalence principles. 




Initial (Cartesian coordinates) New 



Fig. 2. When the material properties are piecewise defined, a pushforward 
of the geometry involving the inverse transformation is useful. 

the cloak). This triangular cylinder is the symmetric of the 
previous one with respect to the horizontal plane containing 
the central fibre of the cylindrical antenna. This bare scatterer 
would therefore give the Fig. 4 image inverted upside-down 
but, here, this object is surrounded by a cloak in order to give 
the very same scattering as before. Indeed, on both sides, the 
electric fields outside the cloak limit are alike. 

Fig. 6 highlights the different scattering patterns by display- 
ing the value of $le(E z ) on a circle of radius 4A located around 
the antenna-scatterer system in the three following cases: the 
case of Fig. 4 (original) with the triangle alone, the case of 
Fig. 5 (coated), and the triangle of Fig. 5 without the coating 
(reversed). It is obvious that the coating restores the field 
distribution independently of the object present in the central 
hole. 

The numerical computation is performed using the finite 
element method (via the free GetDP [16] and Gmsh [17] 
software tools). The mesh is made of 148,000 second order 
triangles including the Perfectly Matched Layers used to 
truncate the computation domain. The singularity of e and [i 
requires a very fine mesh in the vicinity of the inner boundary 
of the cloak (see Fig. 3) and is also responsible for the small 
discrepancies between the numerical model and a perfect cloak 
(see Fig. 5) — including the non zero field in the hole of the 
cloak. 

Note that a small technical problem arises in practice 
when rule (14) is applied: the material properties are defined 
piecewise on various domains and it is very useful to know 
explicitly the boundaries of these domains, e.g. to build 
the finite element mesh (see Fig. 3). These boundaries are 
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Fig. 3. This figure shows a part of the triangular mesh used for the finite 
element modeling of the scattering problem of Fig. 5. The singular behavior 
of the permittivity and of the permeability requires a very fine mesh along 
the inner boundary of the cloak in order to achieve a satisfactory accuracy 
with the numerical model. 




Fig. 4. A conducting triangular cylinder is scattering cylindrical waves. 

curves in the cross section and are thus contravariant objects. 
Therefore, their transformation requires the inverse map (see 
Fig. 2) cp^ 1 from M to TV. Fortunately, map (15) is very simple 
to invert. More explicitly, for a given curve x(f) of parameter 
t in the initial Cartesian coordinates, its push forward by 
Pendry's map is: 

^> = M^ + j^MO, 06) 

with the same variation of the parameter t and ||x|| = r. Note 
that the most common curves used in the design of devices, i.e. 
line segments and arc of circles, are transformed to less usual 
curves except for radial segments (with respect to the center 




Fig. 5. A triangular cylinder different from the one on Fig. 4 is surrounded 
by a cloak designed to reproduce the scattering pattern of the Fig. 4 triangular 
cylinder in spite of the change of scattering object. Of course, the scattering 
object inside the cloak may be arbitrary as far as it is small enough to fit 
inside the cloak. 




reversed 



Fig. 6. The value of the electric field (the real part of E z ) on a circle of 
radius 4A concentric with the cloak is represented as a function of the position 
angle 9 (increasing counterclockwise and with 9 = corresponding to the 
point the most on the right). The three configurations considered here are the 
ones of Fig. 5 (coated), Fig. 4 (original), and the triangle of Fig. 5 without 
the coating (reversed). 

of the cloak) and arc of circles concentric with the cloak. 

On Fig. 5, the image by (p^ 1 of the triangle of Fig. 4 is 
the curvilinear triangle inside the coating region of the cloak. 
In practice, this anamorphosis of the triangle is described by 
three splines interpolating each 40 points that are images of 
points of the segments by ip^ 1 . 

V. Superlens Illusion 

Another dramatic example of transformation optics devices 
are the superlenses [18]: even if these devices were proposed 
a few years before the rise of transformation optics, they are 
nicely interpreted as corresponding to a folding of the space on 
itself. It has been suggested that such devices allow a kind of 
"remote action" of the scatterers making possible things such 
as immaterial waveguides called "invisible tunnels" [19]. They 
can also be used to set up a new kind of invisibility devices 
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[20] and also illusion devices [21] with a similar function that 
the one presented here with generalized cloaking but based on 
negative refraction index materials. This device depends both 
on the object to be transformed, since its scattering pattern has 
to be erased first by an ad hoc "antiobject", and on the object 
to be faked while our device is more general since it depends 
only on the object to be faked independently of the original 
object. As an illustration of superlensing, consider the multiple 
valued transformation of Fig. 7. The part with negative slope 
corresponds to a negative refraction index material (V and 
// have negative eigenvalues) and acts as a superlens. As 
a transformation of empty space, it does not perturbate the 
cylindrical waves emitted by a wire antenna (Fig. 8 where the 
inner disk is the image of four time larger disk) but for the 
attenuation due to the dissipation introduced in the superlens 
permittivity in order to avoid the anomalous resonances [22]. 
One percent of losses has been added here to the value of the 
permittivity of the ideal perfect lens (V has been multiplied 
by 1 — O.Olz). The antenna on the right of the lens has two 
images, one inside the annular superlens and one inside the 
central part of the device so that we have well three copies of 
the antenna. In Fig. 9, a small perfectly conducting deflector 
inside the region surrounded by the perfect lens acts on the 
image of the antenna and forces the waves to propagate only 
to the right. This can also be interpreted as if the deflector has 
a four time larger image acting on the original antenna giving 
the illusion of a much larger object. 




Fig. 7. A superlens can be designed by folding the space on itself, here by 
transforming the radial distance. In this case, there is a part of the physical 
space that has a threefold image in the equivalent space. 



VI. Conclusion 

Transformation optics do not only offer the possibility to 
make optically disappear objects in invisibility cloaks but also 
to completely tune their optical signature i.e. to give them 
an arbitrary appearance. The possibility to place an object 
inside the coating of the cloak and that it will therefore appear 
different was already considered in [3] where a point source 
was shifted creating a mirage effect. A more general case 
is considered here since both the shape and the position of 
the object placed in the coating are modified. Note that if 
the object used to create the illusion is perfectly conducting 
and surrounds the central point of the cloak, its anamorphosis 
will surround the inner boundary of the cloak and it de facto 
suppresses the singular behavior of the material properties. 




Fig. 8. A superlens designed by folding the empty space on itself do 
not perturbate the cylindrical waves emitted by a wire antenna but for the 
attenuation due to the dissipation introduced in the superlens permittivity in 
order to avoid the anomalous resonances [22]. The antenna has an image 
inside the superlens and inside the central part of the device. 




Fig. 9. A small perfectly conducting deflector acts on the image of the 
antenna in the central part of the device and forces the waves to propagate 
only to the right. This can also be interpreted as if the deflector has a four 
time larger image acting on the original antenna. 



In this particular case, the present approach is similar to the 
the hiding under the carpet idea of Li and Pendry [23] but 
here a bounded object replaces the infinite reflecting plane. 
As an illustration of the state of the art cloaking technology, 
this latest device has been recently implemented in the realm 
of infrared using dielectric metamaterials [24]. Note also that 
our technique can be naturally extended to cloaks of arbitrary 
shapes [25], [26]. 

Appendix 
Differential Geometry [5], [6], [9] 

Given a n-dimensional space with a (global) co-ordinate 
system u\ , ■ ■ ■ ,u n (not necessarily orthogonal), the exterior 
derivative d of a function f{u\,---u n ) is its differential 
df = J2i fr^dui. This is a 1-form. A general 1-form can 
be written Y^ i g i (uj)du i where gi(uj) are functions of the 
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coordinates Uj. If a 1-form can be expressed as the differential 
of a function, it is an exact 1-form. 

A curve 7 is an application from an interval I = [t ,ti] 
of K on the n-dimensional space: r(t) = (u\(t), • • • , u n (t)) 
where t is the parameter. The integral fa of a 1-form 
a = gi(uj)dui on tne curve 7 is defined by f^a = 

fto(T,i9i( u j( t ))^^-) dt - The value of the integral depends 
on 7 but does not depend on the choice of the parameter. 

The exterior product A is the skew-symmetric tensor prod- 
uct such that dui A duj = —duj A dui — \{dui ® 
duj — duj ® dui). A general 2-form is a linear combination 
Si j 9ijdui A duj. The exterior derivative of the 1-form 

Y Jl 9idu i is dJ2i9zdu t = J^ij ^ du t A du j- 

A surface S is an application from a two-dimensional open 
domain ft C R 2 on the n-dimensional space: r(s,t) = 
(u\(s, t), ■ ■ ■ ,u n (s,t)) where s,t are the parameters. The 
integral J s f3 of a 2-form /3 = ^ . gijduiAduj on the surface 
£ is defined by 

where d ^" u J> are the Jacobians. The value of the surface 

a(s,t) 

(flux) integral depends on E but does not depend on the way 
the parameters are chosen. 

The Stokes theorem states that J s da = f dJ , a where <9£ is 
the boundary (curve) of the surface S. 

More generally p-forms (with < p < n) are defined as 
totally skew-symmetric tensors and can be manipulated using 
the exterior derivative and the exterior product. Given a 1-form 
a = ctidui + a 2 du 2 + 0:3^3 and a 2-form j3 = fi 23 du 2 A 
du 3 + /?3idu 3 A du\ + fi\ 2 du\ A du 2 , one has for instance: 
da = (^-^)du 1 Adu 2 + (^-^)du 2 Adu 3 + (l^- 
^)du 3 A duu dp = + ^ + |^ )dui A du2 A dtl3> 

and a A f3 = (ax/323 + a 2 $ 31 + a 3 (3 12 )du 1 A du 2 A du 3 . 

All the concepts here above rely only on the topological 
and differential structure of the space. 

The metric is a supplementary structure determined by a 
rank 2 covariant symmetric tensor g whose n 2 coefficients 
form a positive definite matrix. Given a metric, it is possible 
to introduce the concepts of scalar product, norm, distance, and 
angle. The metric allows the definition of a Hodge star operator 
* that is a linear operator on differential forms mapping p- 
forms on (n — p) -forms. 

Particular cases of spaces with a metric are the Euclidean 
spaces E™ where Cartesian coordinates can be chosen so that 
the coefficients of the metric form a unit matrix. For E 3 , 
Cartesian coordinates are denoted {ui = x, u 2 = y, u 3 = z} 
and the metric has the form: g = dx®dx + dy®dy + dz®dz. 
In these Cartesian coordinates, the Hodge operator has the 
following action: 

*dx = dy A dz, *dy — dz A dx, *dz — dx A dy 
*(dx A dy) = dz, *(dz A dx) = dy, *(dy A dz) = dx 
*1 = dx A dy A dz, *(dx A dy A dz) = 1. 
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